Vortex nucleation in Bose-Einstein condensates in an oblate, purely magnetic potential 
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We have investigated the formation of vortices by rotating 
the purely magnetic potential confining a Bose-Einstein con- 
densate. We modified the bias field of an axially symmetric 
TOP trap to create an elliptical potential that rotates in the 
radial plane. This enabled us to study the conditions for vor- 
tex nucleation over a wide range of eccentricities and rotation 
rates. 
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The existence of quantised vortices is one of the most 
striking and fascinating signatures of superfluidity. There 
has been theoretical interest in vortices in a dilute gas 
BEC for many years (for a review see for example Q|) 
and the recent experimental observations |^|| have stim- 
ulated much further work. Vortices have been nucleated 
using phase imprinting (^] and by using a laser beam 
'stirrer' We have used a new, purely magnetic ex- 

citation scheme, which like the optical stirrer has a close 
analogy to the 'rotating bucket' experiment with liquid 
He ||. Vortices provide dramatic visual evidence for 
the single condensate wavefunction, the detailed mech- 
anisms for vortex formation, stabilisation and decay re- 
veal new information about the dynamics of the excited 
states of the condensate and its interaction with the ther- 
mal cloud. The excited states and hence the conditions 
for nucleation depend on the geometry of the confining 
potential |6j. Our work, in an oblate geometry compli- 
ments the work done to date in prolate and spherical 
traps. The range of eccentricity available in our purely 
magnetic rotating trap is much greater than in previ- 
ous experiments, where the radial symmetry was broken 
by 'stirring' the condensate with the dipole force from 
off-resonant laser beams. 

We have measured the relationship between the eccen- 
tricity of the trap in the plane of rotation and the angular 
frequency fl that must be satisfied to nucleate vortices in 
our trap. The trap has a harmonic potential with oscil- 
lator frequencies u> x < ui y < ui z . The shape of the trap is 
characterised by the deformation parameter e, given by 
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We show that when the trap is adiabatically ramped from 
circular, e = 0, to a given eccentricity, there is both an 
upper and lower critical rotation rate for nucleation. We 



present the threshold conditions for vortex nucleation in 
our geometry and finally discuss the role of the thermal 
cloud in both the nucleation and stabilization processes. 

A detailed description of the rotating elliptical trap 
is given in Q] and only a summary of the technique is 
presented here. Our TOP trap consists of a spherical 
quadrupole field and a rapidly rotating bias field The 
zero of the magnetic field describes a circle in the radial 
(XY) plane. The atoms experience a time-averaged, axi- 
ally symmetric harmonic potential. If the X and Y com- 
ponents of the bias field (B x and B y ) are not of equal 
amplitude, then the axial symmetry is broken and the 
potential becomes elliptical in the XY plane. The ra- 
tio of the radial trap frequencies u x /uj y , can be calcu- 
lated numerically from B y /B x . For small eccentricities 
Lo x /Ldy — 1 sa (B y /B x — l)/4. To rotate the elliptical 
potential we modulated B x and B y at 17, a frequency 
much lower than the TOP frequency loq. The final bias 
field has the form of a rotation matrix through angle fit 
multiplying the X and Y components of an ellipse: 
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During any particular evaporative cooling run the value 
of fl, the angular frequency at which the elliptical trap 
rotates, was fixed (for technical reasons) and so the trap 
was always spinning. To create a 'static' trap for efficient 
evaporative cooling, the bias field ratio, E, is set to 1 so 
that the trap is symmetric around the rotation axis. This 
symmetry was checked to be accurate to e = 0±0.005 us- 
ing measurements of the frequency of dipole oscillations 
in the X and Y directions. Evaporative cooling followed 
by an adiabatic expansion gave a condensate of 2 x 10 4 
atoms, at a temperature of 0.5T/T C . At this stage the 
trap is axially symmetric, u> x = oj v , with trap frequencies 
typically ui x /2ir — 62 Hz and w z /2ir = 175 Hz. 

To create vortices, the value of E was ramped linearly 
over 200ms from 1 to its final value, to give a trap that 
was elliptical in the rotating frame. If one creates an ec- 
centric trap by increasing (decreasing) the TOP bias field 
then all three trap frequencies are reduced (increased). 
Since vortex lifetimes depend on the mean trap frequency 
we adjusted the quadrupole field in the adiabatic ex- 
pansion stage, to ensure that all traps had the same av- 
erage radial trap frequency luj_, defined as 
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The condensate was then held in the rotating anisotropic 
trap § for a further 800 ms before being released. After 
12ms of free expansion the cloud was imaged along the 
axis of rotation using an absorption imaging system with 
3 fim resolution. Figure shows images of the expanded 
condensate at different stages during the nucleation pro- 
cess. Initially the cloud elongates, confirming that nu- 
cleation is being mediated by excitation of a quadrupole 
mode (fig.0(a)) [ jlOj ] . Then finger- like structures appear 
on the outside edge of the condensate which eventually 
close round and produce vortices, ~ 800ms after rotation 
began (fig. 0(b)). Approximately 200 ms later, these have 
moved to equilibrium positions within the bulk of the 
condensate. Figures 0(c) and (d) show typical, single- 
shot images of stable vortex arrangements. The depth 
of each vortex (in the integrated absorption profile) is 
approximately 95% of the surrounding condensate. The 
core diameters of ~ 5fim after 12 ms of free expansion 
are consistent with the predictions in The maxi- 

mum number of vortices we observed was seven, limited 
by the number of atoms in our condensate. 

Our first study of the nucleation process involved 
counting the number of vortices as a function of the nor- 
malised trap rotation rate, = Q/u>±, for a fixed ec- 
centricity. Results for trap deformations e = 0.084 and 
0.041 are given in fig. |^These graphs show a maximum 
and minimum value of Q for nucleation at a given eccen- 
tricity. Increasing e increases the range of £1 over which 
vortices may be nucleated, both by lowering Q m in and 
raising VL max . In the limit of small eccentricities, \/2ui± 
is the frequency of the m=2 quadrupole mode, which has 
been shown elsewhere to play a critical role in the nu- 
cleation process Rotation of the ellipsoidal trap at 
half this frequency, i.e. f2 c = 1/V2 ~ 0.71, excites this 
mode. The plots in fig. [^confirm that the nucleation de- 
pends on resonant excitation of the quadrupole mode as 
seen previously |^,|l^l . The resonance is broader at higher 
eccentricity, as intuitively expected for stronger driving. 

Our second study involved holding Q constant and 
counting the number of vortices as a function of the trap 
deformation e. Figure ^ shows our results for two cases: 
(a) f2 > fl c and (b) f2 < 51 c . Interestingly we were 
able to nucleate vortices under adiabatic conditions when 
£1 < £l c . We employ the adiabaticity criterion e/e <C lo±. 
As a further check of the adiabaticity we varied the time 
for the eccentricity ramp between 200 ms and 1 s, and 
detected no difference in the number of vortices formed. 

The critical values of f2 and e for nucleation were ex- 
tracted from plots such as fig. ^| and fig. ^| for many values 
of f2 and compiled on fig. [|. The data points show the 
minimum eccentricity required for a given rotation rate 
and map out region 2, within which vortices nucleate. 
fl c appears to be a critical rotation frequency at which 
vortices can be nucleated with minimum eccentricity, as 



predicted in |12| ] . Changing in either direction requires 
a more elliptical trap for nucleation, although different 
physical processes control the upper and lower limits as 
explained below. 

At f2 > O c , the condensate follows a particular 
quadrupole mode at small eccentricity, region 3 in fig. |^. 
This mode is referred to as the 'overcritical branch' in 
]l3| and has an elliptical density distribution which is 
orthogonal to the trap potential. It then nucleates vor- 
tices when the eccentricity is too large for the quadrupole 
mode to be a solution of the hydrodynamic equations. 
The boundary of the region in the e versus plot where 
this quadrupole mode exists is given by 
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We determined this relation from the solutions of the hy- 
drodynamic equations for superfluids and it is plotted as 
a solid line in fig. || This line agrees well with the exper- 
imental data for the critical conditions for nucleation for 
fl > Q c and a wide range of e. 

Below f2 c , the deformation needed to nucleate vor- 
tices appears to increase linearly with Q. This bound- 
ary cannot be explained in terms of the stability limit 
of a quadrupole mode - the 'normal branch' is stable on 
both regions 1 and 2 and the 'overcritical branch' is sta- 
ble in neither. Our data appears to be at variance with 
the results in |Io| ] , where no vortices were seen when the 
eccentricity was increased adiabatically and fl < fl c . 

A mechanism for the creation of vortices at a frequency 
below f2 c has been proposed in |L2|. They have shown 
that there are regions in the plot of e versus Q, both 
above and below f2 c , where the quadrupole solutions of 
the hydrodynamic equations become dynamically unsta- 
ble. For frequencies above fl c the predicted instability 
domains coincide with the experimentally observed vor- 
tex domains in [fiof , and this work, thus indicating a link 
between their instability analysis and vortices. To make 
a quantitative prediction for the boundary between re- 
gions 1 and 2 shown in fig. ^, will require further detailed 
work for our specific case. We also note that it is possible 
that there exist non-negligible terms of higher order than 
quadratic in the rotating magnetic potential. These have 
been demonstrated to excite high order rotating surface 
modes, e.g. hexapole, which lead to vortex nucleation at 
frequencies less than fi c Jh||| . 

Another possible mechanism for observation of vortices 
below Q c is that the thermal cloud plays an important 
role. Transfer of angular momentum to the condensate 
from the spinning thermal cloud may provide a mecha- 
nism for vortices to form at < S7 C . However the transfer 
rate of angular momentum must be greater than any loss 
rate due to residual trap anisotropy |L4|]. In pQ], grav- 
ity produces a small static eccentricity in the trap in the 
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plane of rotation. The 'spin down' time for a rotating 
thermal cloud in the presence of a deformation parame- 
ter e of only 0.01 is very short, 0.5 s, compared to the 
spin up time of 15 s and hence the cloud may never gain 
significant angular momentum. However, in our experi- 
ment, gravity acts along the rotation axis and hence the 
trap is symmetric in the plane of rotation, giving a more 
favourable ratio of spin-up to spin-down times. 

With this hypothesis in mind, we tested our nucleation 
curve at lower temperature to see if there was any change 
when the amount of thermal cloud was reduced. When 
acquiring the data of fig. [I| the rf was turned off after 
condensation and some heating was observed during the 
nucleation procedure, resulting in a temperature around 
0.8T C . To achieve a lower temperature we left on the so- 
called 'rf shield' so as to give an effective trap depth of 
800 nK during the nucleation process. This resulted in a 
temperature of 0.5T C . No significant change was observed 
in the nucleation curve at this lower temperature. How- 
ever, this does not totally rule out a role for the thermal 
cloud in the nucleation process since even at 0.5T C , there 
was still a significant 20% of atoms in the thermal cloud. 

Although the exact amount of thermal cloud seemed 
to have little effect on the nucleation conditions, it had 
a striking effect on the behaviour of vortices after nucle- 
ation. Without the rf shield during the nucleation pro- 
cess, vortices were only occasionally found in an equilib- 
rium configuration (i.e. 1 vortex in the centre, 3 vortices 
in an equilateral triangle as in fig. [TJ) and - 400ms after 
forming they had already moved to the edge of the con- 
densate before disappearing. However with the rf shield 
present, the vortices were normally found in equilibrium 
positions ~ 200 ms after formation and had a lifetime of 
~ 4s, only limited by the decay of the condensate itself. 

In summary, we have used a purely magnetic rotating 
trap to investigate conditions for vortex nucleation (after 
the formation of the condensate) over a wide range of 
trap eccentricities. For a given eccentricity, we observe 
both an upper and lower limit to the rotation rate for 
nucleation. The upper limit confirms the predictions in 
|iof , but over a much wider range of parameters. However 
the lower limit to the rotation rate is different to that 
reported elsewhere. Further theoretical work is required 
to explain the linear dependence of e on shown in fig. ^. 
The thermal cloud is shown to destabilise vortex arrays, 
thus measurements of the vortex lifetime in an oblate 
geometry were made with an rf shield to prevent heating. 

The rotating anisotropic magnetic trap is well suited 
to experiments that require a large trap eccentricity. It 
has recently been used to investigate the irrotational be- 
haviour of vortex-free condensates at low rotation rates 
|L5j . A large eccentricity is also needed to spin up a ther- 
mal cloud [ pT[ , hence this trap may be useful in attempts 
to condense directly into a vortex state from a spinning 
thermal cloud in thermodynamic equilibrium. 
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(a) 
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FIG. 1. Images of the condensate at different stages dur- 
ing the vortex nucleation process. (All taken with = 0.70 
and e = 0.05 and after 12ms of free expansion), (a) After 
the 200ms spinning eccentricity ramp, the condensate is elon- 
gated, indicating that a quadrupole mode has been excited, 
(b) After a further 600ms in the spinning trap one vortex has 
just formed near the edge. After 800ms, the vortices have 
reached their equilibrium positions and appear in symmetri- 
cal configurations, (c) shows one centred vortex, typical under 
these conditions, whilst (d) shows a triangular array of three 
vortices. 
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FIG. 2. The mean number of vortices as a function of the 
normalised trap rotation rate = tl/u>±. Two different trap 
eccentricities were used, e=0.041 (open circles) and e=0.084 
(solid circles). Each data point is the mean of 4 runs. 



FIG. 3. The mean number of vortices as a function of trap 
deformation at 4 different trap rotation rates: (a) above and 
(b) below the critical value fi c = 0.71. In (a), Q = 0.74 (solid 
circles)and 0.81 (open circles). In (b) SI = 0.61 (solid circles) 
and 0.70 (open circles). Positive (negative) e corresponds to 

LO X < LOy (CJ X > LOy). 
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FIG. 4. The critical conditions for vortex nucleation. The 
data points mark the minimum trap deformation for nucle- 
ation at a particular fi. Vortices may be formed in region 
2. The solid line shows the theoretical limit of stability of 
the quadrupole mode, which is stable in region 3. This line 
is in good agreement with the extreme conditions for vortex 
formation at Q. > fl c 
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